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Abstract. Models for the evolution of dust can be used to predict global 
evolutionary trends of dust abundances with metallicity and examine the relative 
importance of dust production and destruction mechanisms. Using these models, 
we show that the trend of the abundance of polycyclic aromatic hydrocarbons 
(PAHs) with metallicity is the result of the delayed injection of carbon dust that 
formed in low mass asymptotic giant branch (AGB) stars into the interstellar 
medium. The evolution of dust composition with time will have important 
consequences for determining the opacity of galaxies and their reradiated thermal 
infrared (IR) emission. Dust evolution models must therefore be an integral 
part of population synthesis models, providing a self-consistent link between the 
stellar and dust emission components of the spectral energy distribution (SED) 
of galaxies. We also use our dust evolution models to examine the origin of dust 
at redshifts > 6, when only supernovae and their remnants could have been, 
respectively, their sources of production and destruction. Our results show that 
unless an average supernova (or its progenitor) produces between 0.1 and 1 Mq 
of dust, alternative sources will need to be invoked to account for the massive 
amount of dust observed at these redshifts. 



1. Historical Overview 

Chemical evolution (CE) models follow the formation, destruction, abundance, 
and spatial and stellar distribution of elements created during the nucleosynthe- 
sis era of the Big Bang, and the different evolutionary stages in stars. Models are 
pitted against a host of observational test, such as the relative abundances of the 
various elements and their isotopes in meteorites, stars, and in the interstellar, 
intergalactic and intracluster media, the G-dwarf metal hcity distribution, and 



theage- metallicity relation of the various systems [e.g. (jMatteuccil l200ll : iPagel 
19971 )] ■ 



CE models provide a natural framework for studying the evolution of dust 
since the abundance of the elements locked up in dust must be constrained 
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by the availability of refractory elements in the interstellar medium (ISM). CE 
models need then to be generalized to include processes unique to the evolution 
of dust: the condensation efficiency of refractory elements in stellar ejecta, the 
destruction of grains in the ISM by expanding supern o va rem nants (SNRs), and 
the growth and coagulation of grains in clouds iDw eg (199 ^). 

The first dust evolution models (|Dwek fc Sc alo 1979 . OdSdi ) (hereafter DS) 
addressed the origin of the eleme ntal depletion pattern, which was a subject of 
considerable debate. Field ( 19741 ) showed that the depletion pattern correlated 
with condensation temperature, suggesting that it reflects the condensation ef- 
ficiency of the elements in their respective sources. Such causal correlation 
requires that dust undergoes very little interstellar processing that can alter the 
depletion pattern. An equally good correlation exists bet ween the ma gnitude of 
the elements' depletion and their first ionization potential ( Snow|[l975l ) . suggest- 
ing that the depletion pattern may instead be governed by accretion processes 
in molecular clouds. Pointing out the intrinsic physical correlation between the 
condensation temperature, the first ionization potential, and the threshold for 
grain destruction by sputtering, DS suggested that the depletion pattern could 
refiect the destruction efficiency of dust in the ISM. It is currently clear that 
all three processes play some role in establishing the elemental depletion pat- 
tern, since it depends on the density of the medium in which it is observed 
( Savage &i Sembachlll996l ). Globally, their relative importance depends on the 
prevalence of and the cycling times between the different phases (hot, neutral, 
and molecular) of the ISM. 

A more detailed review of dust evolution models is given bv lDweld (1 19981 ). 

Since then, sever al models have been const ructed to follow the evol ution of dust 

in dw arf galaxies ( Lisenfeld Sz Ferraralll998l). Damped Lya systems ( Kasimova Sz Shchekinov 
20031 ). the Milky Way ga laxy (I Zhukovska et al."2008'), and the origin of isotopic 



anomalies in meteorites (jClavton &: Nittler 2004 : .Zinner et al.i . 2006a, b). 

Early dust evolution models adopted the instantaneous recycling approxi- 
mation, which assumes that all the elements and dust are promptly injected back 
i nto t he IS M following the formation of their nascent stars. In contrast iDweS 
( 19981 ) and [Morgan Sz Edmundsl ( 20031 ) constructed models that take the finite 
main sequence (MS) lifetimes of the stars into account. Silicate dust is primarily 
produced in supernovae (SN) that "instantaneously" recycle their products back 
to the ISM, whereas carbon dust is mainly produced in low-mass carbon-rich 
stars which have significantly longer MS lifetimes. Consequently these models 
predicted that the composition of the dust should evolve as a function of time. 

Dust evolution models contain many still uncertain parameters such as the 
dust yields in the various sources and the grain destruction efficiency in the 
ISM. However, in spite of these uncertainties we will show that they can be very 
successful in predicting global evolutionary trends, namely the observed correla- 
tion of the abundance of polycyclic aromatic hydrocarbon (PAH) molecules with 
metallicity, and in examining the origin of dust in the high-redshift universe. 
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2. Main Ingredients of Dust Evolution Models 
2.1. The Yield of Dust in Stars 

The dust condensation efficiency, and its composition and size distribution de- 
pend very much on the environment in which it is formed. Dust can form in the 
quiescent outflows of AGB or Wolf-Rayet stars, or in the explosively expelled 
ejecta of SNe and novae. There is substantial observational evidence for the for- 
mation of dust in all these sources, however, their relative importance as sources 
of interstellar dust, and the composition and size distribution of the condensed 
dust are still uncertain, especially in SNe. Because of the ease of formation and 
stability of the CO molecule, dust sources can be divided into two categories: 
carbon-rich sources, in which the C/0 abundance ratio is larger than 1. These 
sources will produce carbon dust; and oxygen-rich sources with C/0 < 1, which 
will produce silicate-type dust. Such simple arguments assume that CO forma- 
tion goes to completion, which ma y not always be the case , as suggested by the 
dynamical condensation models of iFerrarotti &: Gail ( 20061 ). 



SN explosions mark the death of stars more massive th an ^ 8 M(^. Fi gure 
[U depicts the post-explosive composition of a 25 Mq star (IWooslevI Il988l ). It 



depicts a typical onion-skin structure in which the composition of the differ- 
ent layers reflect the pre- and post-explosion nuclear burning stages of the star. 
Globally, the ejecta has a C/0 ratio < 1, and should in principle be only produc- 
ing silicate dust. However, in spite of the mixing between the layers caused by 
Ray leigh- Taylor instabilities in the ejecta, this mixing is of macroscopic nature 
and does not occur on an atomic level. So the layers above above ~ 4.2 Mq, in 
which C/0 > 1, will maintain this ratio and produce carbon dust, whereas the 
inner layers will produce silicates. If all condensible elements precipitated out 
of the ejecta and formed dust, t he yield of dust in a typical SN would be about 
1 Mq (|Kozasa et al.lll989l . Il99lh . 



The yield and composition of dust in lower mass stars depends on the 
C/0 ratio in their atmosphere during the AGB phase of their evolution. Fig- 
ure [2] depicts the C and O yield in stars. Supernovae yie lds were taken from 



Wooslev k Weaver! (|l995l ) and AGB yields were taken from feara kas &: Lattanziol 



20031 ). Stars with masses above ~ 8 Mq produce carbon and silicate dust, irre- 



spective of the global C/0 ratio in their ejecta. The mass range of carbon rich 
stars depends on the initial stellar metallicity. At zero metallicity (left panel) 
the mass range of stars producing carbon dust is between 0.8 and 7.8 Mq. This 
range narrows significantly to masses between 3.0 and 3.6 Mq at solar metal- 
licity (right panel). Figure [3] presents a qualitative depiction of the method we 
use to calculate the yield of carbon and silicate dust in AGB stars. A more 
realistic model calculat ing the yield of dust in AGB stars was presented by 
Ferrarotti Gaill (|2006l ^ and Hoefner (2009, this confer ence proceedings). 



Figure H] (left panel) depicts the stellar evolutionary tracks of stars with an 
initial solar metallicity. Also shown in the figure is t heir MS lifetime. Sim ilar 
lifetimes are obtained for different initial metallicities ( Portinari et al. 19981 ). At 



metallicity of Z = 0, the first carbon dust producing stars are about 8 Mq and 
will evolve of the MS about 50 Myr after their formation. At solar metallicities, 
the production of carbon by AGB stars will be delayed by about 500 Myr, when 
~ 4 Mq stars evolve off the MS. Since most of the interstellar carbon dust is 
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Figure 1 . The post-explosive composition of a 25 Mq star (|Woosleylll988| ) . 
A typical SN can potentially produce 1 Mq of dust. 




Figure 2. The C and O yields in stars. Stars with masses above ^ 8 Mq 
become SNe and produce both silicate and carbon dust. Stars with masses 
below ~ 8 Mq produce either carbon or silicate dust, depending on the C/0 
ratio in their ejecta. The light gray area in the horizontal bar depicts the 
range of stellar masses in which C/0 > 1 . 



Carbon 
dust 



Silicate 
dust 




Mg Si Fe 



Mg 



Fe 



Figure 3. Qualitative depiction of the calculated yield of carbon and silicate 
dust in AGB stars. When C/O > 1 (left panel), the star produces only carbon 
dust. The dark shaded area depicts the number of carbon atoms that condense 
into dust. When C/O < 1 (right panel) the star produces silicate dust, and 
the dark shaded area depicts the condensing elements. 



made in AGB stars, this delay in its formation can have important observational 
consequences. Carbon dust has a significantly higher visual opacity than sili- 
cates, so the opacity of galaxies will change with time, with young systems being 
more transparent than older ones. 

Figure H] (right panel) depicts the different evolutionary trends of SN- and 
AGB-condensed dust calculated for a CE model with by exponential star for- 
mation rate characterized by a decay time of 6 Gyr, and a Salpeter initial mass 
function (Dwck 1998.. .2005i ). The silicate and carbon dust yields were calculated 
assuming a condensation efficiency of unity in the ejecta, and grain destruction 
was neglected. The model therefore represents an idealized case, in which grain 
production is maximized, and grain destruction processes are totally ignored. 
Also shown in the figure are the separate contributions of AGB stars to the 
abundance of silicate and carbon dust. The onset of the AGB contribution to 
the silicate abundance starts at i ~ 50 Myr, when ~ 8 Mq stars evolve off the 
main sequence, whereas AGB stars start to contribute to the carbon abundance 
only at t ~ 500 Myr, when 4 Mq stars reach the AGB phase. The figure also 
presents the dust-to-ISM metallicity ratio, which is almost constant at a value 
of ~ 0.36. 



3. The Lifetime of Interstellar Dust 



Following their injection into the ISM, the newly-formed dust particles are sub- 
jected to a variety of interstellar processes resulting in the exchange of elements 
between the solid and gaseous phases of the ISM, including: (a) thermal sput- 
tering in high-velocity (>200 km s~^) shocks; (b) evaporation and shattering by 
grain-grain collisions in lower velocity shocks; and (3) accretion in dense molec- 
ular clouds. Detailed description of the various grain destruction mechanisrn s 
and grain lifetimes in the ISM were presented by ( Jones et al. 19961 : Joned 2003 ) . 
In addition, SN condensates can be destroyed shortly after their formation by 
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Figure 4. Left panel: The H-R diagram of stars and their main sequence 
hfetime. Right panel: The evolution of sihcate (dashed line) and carbon (solid 
line) dust from SNe (bold curves) and AGB stars (light curves). 



reverse shocks that travel through the expanding ejecta (Dwek 2005, Bianchi &: 
Schneider 2007, Nozawa et al. 2008). 

The most important parameter governing the evolution of the dust is its 
lifetime, t^, against destruction by SNRs. In an interste llar medium with a 
uniform dust -to-gas mass ratio, Z^^, this lifetime is given by (jPwek fc Scaldll980l : 
McKeelll989l ): 

Mrf(t) _ Mg{t) 



Td 



{rrid) RsN {misu) RsN 



where and Mg are, respectively, the total mass of dust and gas in the galaxy, 
(rrid) is the total mass of elements that are locked up in dust and returned by 
a single SNR back to the gas phase either by thermal sputtering or evaporative 
grain-grain collisions. Rsn is the SN rate in the galaxy, so that the product 
(m^) Rsn is the destruction rate of dust in the ISM. The parameter (misM) = 
(md) jZd is the effective ISM mas s that is completely cleared of dust by a single 
SNR, given by (|Dwek et al.ll2007l ): 



("iISM) 



Do 



dM 



dv. 



(2) 



where Cdi^s) is the fraction of the mass of dust that is destroyed in an encounter 
with a shock wave with a velocity Vs, {dM / dvs)dvs is the ISM mass that is swept 
up by shocks in the [vs, Vg + dvg] velocity range, and and Vf are the initial 
and final velocities of the SNR. Figure [5] depicts the mass fraction of carbon 
and silicate dust that is destroyed after being swept up by a shock of velocity 
Vs as a function of shock velocity. An updated version for the carbon and 
PAH destruction efficiency was presented by Jones et al. (2009, this conference 
proceedings). For example, in the Milky Way Mg 5 x 10^ Mq, Rsn ~ 
0.03 yr~^, and (misM) ~ 300 M©, giving a dust lifetime of ~ 6 x 10*^ yr. 
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Figure 5. The mass fraction of carbon and silicate dust that is destroyed 
as a function of shock velocity (after Jones et al. 1996). 



In addition, SN condensates can be destroyed shortly after their formation 
by reverse shocks that travel through their expanding ejecta (Dwek 2004, Nozawa 
et al. 2008). 



4. PAHs and Silicate Dust as Tracers of AGB- 
Dust 



and SN-condensed 



An exciting discovery made by spectral and photometric observations of nearby 
galaxies with the Infrared Space Observatory (ISO) and Spitzer satellites was 
the striking correlation between the strength of their mid-infrared (IR) aromatic 
features, commonly attributed to the emission from PAHs, and their metallicity, 
depicted in Figure E] [left panel; see iGalliano et all ( 20081 ) for references]. The 
figure shows the rise of i^8/24) the 8 //m-to-24/im band flux ratio with galaxies' 
metallicity, and the existence of a metallicity threshold below which ^3/24 is 
equal the flux ratio of the dust continuum emission. The strength of the aromatic 
feature is a measure of the PAH abundance. Since PAHs are predominantly made 
in C-rich AGB stars, this correlation provides the first observational evidence 
for the delayed injection of AGB condensed dust into the ISM, provided the 
metallicity is a measure of the galaxies' age. The testing of this hypothesis 
requires first the determination of the PAH abundance in each galaxy. 

PAHs are very small macromolecules, typically 50 A in diameter, that are 
stochastically heated by the ambient radiation field. Consequently, only a frac- 
tion of the PAHs are radiating at mid-IR wavelengths at any given time. To 
determine the total abundance of PAHs, including those too cold to emit in the 
aromatic features requires the determination of the intensity of the interstellar 
radiation field (ISRF) to which they are subjected. 
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Figure 6. Left panel: The observed correlation between the 8-to-24 /im 
bands flux density and metaUicity. Right panel: PAH and dust abundances 
derived from detailed models of galaxies' SED versus galaxies' metallicity. 



Figure [7] depicts the steps used by iGalliano et al. I (120081 ) in modeling the 
galaxies' spectral energy distribution (SED). The galaxy used for this illustrative 
purpose is the starbu rst M82. The dust model used in the calculations is the 
BARE-GR-S model of lZubko et aD (|2004l ). consisting of PAHs, and bare silicate 
and graphite grains with solar abundances constraints. 

Figure 7a shows the SED of M82, and its various emission components: 
stellar emission at optical, and near-IR wavelengths; the PAH spectrum at mid 
IR wavelengths; dust emission from mid- to far-IR wavelengths; and free-free 
and synchrotron emission at radio wavelengths. 

Fig 7b depicts a fit to the dust spectrum using an ISRF characterized by a 
power-law distribution of radiation field intensities. PAH abundance determined 
by this method will underestimate the real abundance of PAHs in the galaxy 
compared to the method outlined below. Our models use the free-free and mid- 
IR emissions to constrain the gas and dust radiation from the gas and dust 
from H H regions, and the far-IR and optical emission to constrain the ISRF 
that heats the dust in photo-dissociation regions (PDRs). The radio emission is 
uniquely decomposed into free-free (dashed) and synchrotron (dotted) emission 
components (see Figure 7d). 

Massive stars are required to produce the ionizing radiation and the ex- 
panding SN blast waves that generate, respectively, the observed free-free and 
synchrotron emission. These stars are produced in an "instantaneous" burst of 
star formation, in contrast to the stars that are continuously created over the 
lifetime of the galaxy and mostly contribute to the optical and near-IR emission 
(Figure 7c). 

The ionizing and a fraction of the non-ionizing radiation emitted by the 
starburst component that is absorbed in the H II region is shown as a shaded 
area in Figure 7d. This energy is r eradiated by the dust and gas, giving rise to 
the thermal IR and free-free emission components shown in the figure. 
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The non-ionizing radiation from the older stellar population and the radi- 
ation escaping the H II regions form the diffuse ISRF that is absorbed by the 
dust in photodissociation regions (PDR) (Figure 7e). The shaded area depicts 
the fraction of the radiation from the older stellar population that is absorbed 
in PDRs. The absorbed radiation is reemitted by the dust, giving rise to the IR 
emission (figure 7f). 

Figure 7g depicts the different emission components, and Figure 7h shows 
the fit of th eir sum to the SEP of M82. The details of the fitting procedure are 
described in lGalliano et all ( 20081 ). 

Using this physical fitting procedure, iGalliano et al.l ( 20081 ) derived the 
abundance of PAHs, silicates, and graphite grains in 35 nearby galaxies with 
metallicities ranging from 1/50 to 3 times solar. The detailed physical modeling 
of their SEDs gives larger PAH abundances compared to models that employ a 
template ISRFs to heat the PAHs and the dust. In these models, such as the 
one depicted in Fig. 7b, PAHs are subjected to the same intense radiation field 
as that required to produce the mid-IR emission from hot dust. In our model, 
PAHs are subjected to a weaker radiation field. Since PAHs do not survive 
in H II regions, all their emission originates from PDRs, which are subjected 
to lower intensity radiation fields than the H II regions. So, compared to the 
template ISRF models, a larger amount of PAHs is required to produce the 
same PAH spectrum with a weaker ISRF. The resulting PAH abundances are 
plotted versus metallicity in Figure [6] (right panel). The figure shows that the 
observed trend of increasing 8-to-24 //m band flux ratio with metallicity indeed 
reflects a trend of increasing PAH abundance with metallicity. The figure also 
shows the distinct evolutionary trends of PAHs and the far-IR emitting dust 
with metallicity. 

Figure [8] compares the evolution of the PAHs and dust components derived 
from the dust evolution model to the trend of PAH abundances with metallicity. 
For the sake of this comparison, the evolution of dust abundances as a function 
of time was converted to an evolution as a function of metallicity, using the age- 
metallicity relation derived in the model. We emphasize that the parameters 
used in the dust evolution model (the star formation rate, the stellar IMF) were 
identical to those used in the population synthesis model that was used to fit 
the galaxies' SED. From the dust evolution model we already derived the two 
distinct evolutionary trends of SN- and AGB-condensed dust (see Fig. O right 
panel for the idealized example). The current figure compares these results with 
the derived abundance of the PAHs and of the dust that gives rise to the far- 
IR emission. The latter is dominated by emission from silicates, and should 
therefore follow the trend of the SN condensates, since most silicate grains are 
produced in SNe. The figure shows that the observed far-IR emitting dust falls 
on the evolutionary track of the calculated SN-condensed dust, and that the 
observed PAH abundances fall on the evolutionary track of the carbon dust that 
formed in AGB stars. The shaded regions in the figure represent the range of 
evolutionary tracks that correspond to different parameters that determine the 
star formation rate and grain destruction efficiencies in the models. 



10 




11 




5. The presence of massive amounts of dust at high redshift 

The detection of massive amounts of dust in hyperluminous IR galaxies at 
redshifts z > 6 raises challenging questions about the sources capable of pro- 
ducing s uch large am ount of dus t during the relat ively short lifetime of these 
galaxies (jMaiolino et a l. 2006; B eelen et a l. 2006; iMorgan fc Edmund s 2003). 



For example, the galaxy SDSS J1148+5251 (hereafter J1148+52 51) located at 



z = 6.4 was observed at far-IR and submi llimeter wavelength (jBertoldi et al.l 



20031 : iRobson et alll2004l : iBeelen et al1l2006l ). The average IR luminosity of the 



source is Ltr ^ 2 x 10^^ L tD, and the average dust mass is Md ~ 2 x 10^ Mq. 
Using the Kennicutt ( 19981 ) relation, one can derive a star formation rate (SFR) 



of ~ 3000 Mq yr~^ from the observed far-IR luminosity. For comparison, the 
Milky Way galaxy is about 10 Gyr old, has an average SFR of ~ 3 Mq yr^^, and 
contains about 5 x 10^ Mq of dust, a significant fraction of which was produced 
in AGE stars. 

At z = 6.4 the universe was only 890 Myr old, using standard ACDM 
parameters (Q^ = 0.27, Qa = 0.73, and Hq = 70 km Mpc~i). If J1148-h5251 
formed at z = 10 then the galaxy is only 400 Myr old. If the SFR had occured 
at a constant rate over the lifetime of the galaxy, its initial mass should have 
been about 10^^ Mq, which is si gnificantly larger t han the dynamical mass 
Mdyn ~ 5 X 10^° Mq of the galaxy (jWalter et aDl2004l 1. The high observed SFR 



may therefore represent a recent burst of star formation that has lasted for only 
about 20 Myr. The galaxy J1148-F5251 is therefore at most ~ 400 Myr old, 
and probably significantly younger with an age of only ~ 20 Myr. Adopting 
a current gas mass of Mg = 3 x lO^'^ Mq for this galaxy we get that the gas 
mass fraction at 400 Myr is about 0.60. The dust-to-gas mass ratio is given by 
Zd = Md/Mg = 0.0067. 

A significant fraction of the dust in the Milky Way was produced in AGB 
stars. However, these stars are not likely to contribute significantly to the for- 
mation of dust in very young galaxies, since the low mass stars (M ~ 3 Mq) 
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Figure 8. Comparison between the metallicity trends of the PAH abun- 
dance derived from the observed SED and those derived from the chemical 
evolution model. The shaded area represent the range of model prediction for 
different gra in destruction and sta r formation rates. Details of the figure are 
described in iGalliano et al.l llOOl). 
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that pr oduce most of the dust did not have time to evolv e off the main se- 
quence (iDwekl 119981 : iMorgan fc Edmunds! [200^ : iDwekl l2005l ^ . In contrast, core 
collapse SNe (M > 8 Mq) and their post-main-sequence progenitors inject their 
nucleosynthetic products back into the ISM shortly (t < 20 Myr) after their 
formation, resulting in the rapid enrichment of the interstellar medium (ISM) 
with the dust that formed during the mass loss phase prior to the SN event, or 
in the explosive SN ejecta. We will hereafter attribute both contribution to the 
SN event, since both are return "promptly" to the ISM. But can SNe account 
for the large amount of dust seen in this object? The answer to this question is 
complicated by the fact that SNe are also t he main source of grain destr uction 
during the remnant phase of their evolution ( Jones et al. 19961 : J ones! 2003 ) . The 
problem can therefore only be quantitatively addressed with CE models for the 
dust in these systems. 

The results of the detailed dust evolution models described summarized in 
Figure [8] show that the contribution of AGB stars to metal and dust abundance 
can be neglected in galaxies with ages less than about 400 Myr. The equations 
for the chemical evolution of the galaxy can then be considerably simplified using 
the instantaneous recycling approximation, which assumes that stars return their 
ejecta back to the ISM promply after their formation . The evolution of the dust 
abundance can t hen be written in analytical form ( Morgan &: Edmundsl 20031 : 
Dwek et al.l[2007l l. 

In particular, the yield of dust, Y^, required to obtain a given dust-to- gas 
ma ss ratio, Z^, when the galaxy reaches a given gas mass fraction ^g, is given 



by (iDwek et al.ll2007l ): 



{msM) 



R 



1 



(3) 



where, 



(1 - R) 



(4) 



R is the fraction of the stellar mass that is returned back to the ISM during the 
stellar lifetime, (misM) is given by eq. (2), and is the mass of all stars born 
per SN event. For example, = 147 and 50 Mq, respectively, for a Salpeter 
and top-heavy IMF. 

Figure [9] shows how much dust an average SN must produce in order to give 
rise to a given dust-to-gas mass ratio, for various grain destruction efficiencies. 
The value of Yd was calculated when jjLg reaches a value of 0.60, the adopted 
gas mass fraction of Jl 148+5251 at 400 Myr. Calculations were performed 
for two different functional forms of the stellar IMF: a Salpeter IMF in which 
(f){rn) ~ m~-^'^^ and 0.1 < m{MQ) < 100; and a top heavy IMF characterized 
by the same mass limits but a flatter slope (/)(m) ~ m~^'^^. Here, (/){m) is the 
number of stars per unit mass interval, normalized to unity between 0.1 and 

100 Mq. 

The figure shows that, for example, to produce a value of = 0.0067 
at fig= 0.60, a SN must produce about 0.4 (1.2) Mq of dust for a top-heavy 
(Salpeter) IMF, provided the dust is not destroyed in the ISM, that is, (th-ism) 
= 0. Even with modest amount of grain destruction, (mispvi) = 100 Mq, the 
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Figure 9. The IMF-averaged yield of dust by type II supernova, Yd, that 
is required to account for a given dust-to-gas mass ratio Zd-, is presented for 
different values of (toism) given in units of A/q. Solid and dashed lines corre- 
spond to calculations done for a top-heavy and a Salpeter IMF, respectively. 
The horizontal dashed line near the bottom of the figure corresponds a value 
of Y d = 0.054 M(n, th e highest inferred yield of dust in supernova ejecta to 
date lRho et al.l ()2008D . The vertical dotted line represents the value of Zd at 
[Xg = 0.60. Curves are labeled by (misM) given in units of Mq. The top two 
dashed (solid) horizontal lines represent IMF-averaged theoretical dust yields 
for a Salpeter (top-heavy) IMF, assuming 100% condensation efficiency in the 
SN ejecta. 



required SN dust yield is dramatically increased to about 1 — 2 M©, depending 
on the IMF. The horizontal line in the figure corresponds to a value of 1^ = 
0.054 Mq, the large st mass of SN-c o ndensed dust inferred to be present in a 
supernov a or SNR (iRho et al.1 l2008l : ISueerman etaP l2006l ). Contrary to the 
claim by iRho et aL ( 20081 ). this yield is not sufficient to account for the large 
amount of dust ob served in high redshi ft galaxies, since the quoted chemical 
evolution models of Morgan Edmundsl ([2003) do not include the effect of grain 
destruction. The figure shows that even without grain destruction, the largest 
observed yield can only give rise to a dust-to-gas mass ratio of ~ 4 x 10~^. If 
the mass of dust in the ejecta of Gas A represents a typical SN yield, then other 
processes, such as accretion onto preexisting grains in molecular clouds is needed 
to produce the mass of dust in J1148-I-5251. 



15 



6. Summary 

Dust evolution models have proven to be very successful in predicting global 
evolutionary trends in dust abundance and composition, and in analyzing the 
origin of dust in the early universe. An important prediction of these models 
is that SN- and AGB-condensed dust should follow distinct evolutionary paths 
because of the different stellar evolutionary tracks of their progenitor stars. By 
analyzing the UV-to-radio SED of 35 nearby galaxies we have identified silicates 
and PAHs, respectively, as tracers of SN- and AGB-condensed dust. Our SED 
fitting procedure used chemical evolution, dust evolution, and population syn- 
thesis models in a consistent fashion. The models used the free-free and mid-IR 
emissions to constrain the gas and dust radiation from the gas and dust from 
H II regions, and the far-IR and optical emission to constrain the ISRF that 
heats the dust in PDRs. The observed correlation of the intensity of the mid-IR 
emission from PAHs with their metallicity can then be interpreted as the result 
of stellar evolutionary effects which cause the delayed injection of carbon dust 
into the interstellar medium. 

The early universe is a unique environment for studying the role of massive 
stars in the formation and destruction of dust. The equations describing their 
chemical evolution can be greatly simplified by using the instantaneous recycling 
approximation, and by neglecting the delayed contribution of low mass stars to 
the metal and dust abundance of the ISM. Neglecting any accretion of metals 
onto pre-existing dust in the interstellar medium, the evolution of the dust is 
then primarily determined by the condensation efficiency of refractory elements 
in the ejecta of Type II supernovae, and the destruction efficiency of dust by SN 
blast waves. 

We applied our general results to J1148-I-5251, a dusty, hyperluminous 
quasar at redshift z = 6.4 and found that the formation of a dust mass fraction 
of Zd = 0.0067 in a galaxy with an ISM mass of 3 x lO^'' Mq, requires an average 
SN to produce between 0.5 and 1 Mq of dust if there was no grain destruction. 
Such large amount of dust can be produced if if the condensation efficiency 
in SNe is about unity. Observationally, the required dust yield is in excess of 
the largest amount of dust (~ 0.054 Mq) observed so far to have formed in a 
SN. This suggests that accretion in the ISM may play an important role in the 
growth of dust mass. For this process to be effective, SNRs must significantly 
increase, presumably by non-evaporative grain-grain colUsions during the late 
stages of their evolution, the number of nucleation centers onto which refractory 
elements can condense in molecular clouds. 
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